Nonlinear anomalous diffusion equation and fractal dimension: Exact generalized 

gaussian solution 
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In this work we incorporate, in a unified way, two anomalous behaviors, the power law and 
stretched exponential ones, by considering the radial dependence of the TV- dimensional nonlinear 
diffusion equation dp/dt = V ■ (KV p v ) — V ■ {pFp) — ap, where K — Dr~ 6 ', v, 8, p and D are real 
parameters and a is a time-dependent source. This equation unifies the O'Shaugnessy-Procaccia 
anomalous diffusion equation on fractals (y = f) and the spherical anomalous diffusion for porous 
media (0 — 0). An exact spherical symmetric solution of this nonlinear Fokker-Planck equation 
is obtained, leading to a large class of anomalous behaviors. Stationary solutions for this Fokker- 
Planck-like equation are also discussed by introducing an effective potential. 

PACS numbers: 05.20.-y, 05.40.Fb, 05.40.Jc 
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I. INTRODUCTION 

Anomalous diffusion has been the subject of numerous 
investigations in the last years. In particular, anomalous 
diffusion has played a fundamental role in the analysis of 
a wide class of systems such as diffusion in plasmas [jj], 
turbulent flows 0, transport of fluid in porous media 
p|, and diffusion on fractals H. In general, anomalous 



J 



diffusion may be classified employing the mean square 
displacement, ((Aa:) 2 }. When {(Ax) 2 ) oc t v , the value 
X] > 1 characterizes a superdiffusive process, rj < 1 a 
subdiffusive one and 77 = 1 a 'normal' diffusion. 

We consider here a quite large class of anomalous dif- 
fusion, namely those associated with the ./V-dimcnsional 
nonlinear equation Q 
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dp(r, t) 
dt 



= V • {KS/[ P (r, t)]"} - V • [pF(r, t)p(r, t)} - a(t)p(r, t), 



(1) 



where K = Dr~ , v, 6, and D are real parameters, p, 
is the drift coefficient, F(r,t) is the external force, and 
a(t) is a time-dependent source. From the Gauss theorem 
one verifies that this general equation preserves the norm 
J p{r, t)d N r if a(t) = and KV[p{r, t)]» - pF(r, t)p{r, t) 
goes sufficiently rapid to zero when r = |x| — > 00. The 
difference between the Fokker-Planck-like equation (|l|) 
and the usual one is that the diffusion term depends 
on a power of p(r, t) and the diffusion coefficient is spa- 
tially dependent. Notice that this equation for = and 
a(t) = has the form: 



dp(r,t) 
dt 



Z)V 2 [p(r,^-V-[ M F(r,t)p(r,i)], (2) 



and it is usually called porous media equation. The 
above equation admits important physical applications 
such as percolation of gases through porous media 
thin liquid films spreading under gravity 0, some self- 
organizing phenomena and surface grow [^|, among 
others (see Ref. || and references therein) . Furthermore, 
this equation has been related to nonextensive Tsallis 
statistics. This connection was first pointed by Plas- 



tino and Plastin o |l0| a nd revisited by others authors 

Another important kind of anomalous diffusion 
emerges from Eq. ([!]) when v = 1, a(t) = and 
F(r, t) = . In this case, the equation 



dp{r, t) 
dt 



V-[KV P (v,t)\ 



(3) 



was related to turbulent diffusion in the atmosphere 
@ (see also Ref. @). Here K oc r 4 / 3 . In gen- 
eral, in a (i-dimensional space, V • K\J is proportional 



to 



-1-6 d 



+ 



A 



1 dr dr 

pends on the angular variables). Thus, 



(A is an operator, which de- 



A = 



1 d 



dr 



(4) 



is the radial part of the operator V • KV to be considered 
in the study of d-dimensional spherical symmetrical solu- 
tions of Eq. (|3|) . If d is interpreted as a positive real num- 
ber, it plays a role of a fractal dimension embedding in a 
A^-dimensional space. Thus, Eq. (||) becomes, in a fractal 
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framework, the O'Shaughnessy-Procaccia equation) 19 
dp(r,t) 



at 



= DAp(r,t). 



(5) 



We are ready now to unify Eqs. (Q) and (||). In fact, 
if we consider only diffusive terms and the radial part of 
these equations we have [pp| 



dp(r, t) 

at 



= DA[p(r,t)y. 



(6) 



Of course, if = and v = 1, the usual diffusion 
equation is recovered, i.e., 



dp(r,t) 
dt 



DV 2 P (r,t). 



(7) 



In the previous scenario, Eq. (|l|) can be generalized 
in order to be also employed in a fractal context, when 
spherical symmetric solutions are considered and V • K\J 
is replaced by A. This paper is dedicated to investigate 
a class of these solutions in the presence of a source term 
and with a drift term, as well as the stationary solutions 
when a general drift term is present. 

The present work is organized as follows. In the next 
section we review an exact solution (generalized gaus- 
sian) for Eq. (^|) and its limit cases. This generalized 
gaussian essentially gives the basis for the remaining so- 
lutions investigated here. In Sec. Ill, a formal solution 
for the nonlinear Fokker-Planck equation Eq. (||) with a 
source term is obtained. Sec. IV deals with the applica- 
tion of our formal solution for a time power law source. 
In Sec. V, we present a generalized gaussian solution for 
the nonlinear Fokker-Planck with a linear spherical force 
term. Furthermore, in Sec. VI, we consider the station- 
ary solutions for nonlinear Fokker-Planck equation with a 
general drift term. A summary of our results is presented 
in the last section. 



II. UNIFIED POWER-LAW AND STRETCHED 
EXPONENTIAL SOLUTION 

Motivated by the Gaussian solution of the usual diffu- 
sion, Eq. (Q), it is natural to employ a generalization of 
such solution when we are investigating Eq. (|Tj|). In this 
direction, we follow Ref. [pFjj employing the ansatz: 



where A = DX{2-q)[X + d(l-q)}(3 Z^\ (3 = (3(0) and 
Zq = Z(Q). It is also important to note that the product 
(3Z x l d is time independent. 

From the solution (|J) the mean value for a generic 
function f(r) can be calculated. In particular, the mean 
value for r° is given by 



JT V p(r,t)r d - 1 dr 
Jo°° P( r > t)r d ~ l dr 



with 



G a = 



£(fg) 
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(?-ir /A r(^- 
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q<l 
q>l. 



(11) 



(12) 



Observe that the existence of (r a ) imposes restrictions 
over the parameters: A + d(l — q) > a(q — 1) and A > 0. 
In addition, the restriction q < 2 (y > 0) is necessary for 
p(r, t) to be real. In the following, we assume that the 
parameters obey the above restrictions. 

An important particular case of (r a ) is the mean square 
displacement, (r 2 ), since such mean value is commonly 
employed to classify the diffusion process. Thus, by tak- 
ing a = 2 in Eqs. ( |lT| ) and (|l^) we have 



(r 2 ) = C 2 f3- 2/X {l + At] 



i2/[A+d(l- g )] 



(13) 



Consequently, (r 2 ) ~ t 2 /l' 2 + e + d ^- 1 )} for a sufficient large 
t. Now, it is easy to verify that 9 = d(l — v) describes 
a 'normal' diffusion, even when 9 ^ and v ^ 1. In 
fact, there is a competition between 9 and v in such way 
that the anomalous diffusive regime induced by 9 ^ is 
compensated by a convenient one with v ^ 1. Further- 
more, this competition leads to a subdiffusive process 
when 9 > d(l — v), and to a superdiffusive one when 

e < d(i - v). 

Note, also, that an alternative point of view to verify 
how much a solution deviates from the usual one is to 
analyze the decay of p(0,t). In other words, the way 
the solution decreases in this central point can inform 
about its deviation from the normal diffusion. In our 
case, Eq. (§) gives p(0, t) - t -d/[i+e+d{v-X)} in contrast 
with p(0,t) ~ £~ d / 2 , obtained from the usual diffusion. 

We remark that important limiting cases can be ob- 
tained from Eq. (@). If 9 = 0, the solution for the porous 
media equation without the drift term is recovered j^, [TJJ , 



P(r, t) = [1 - (1 - q)(3(t)r X } 1/(1 q) I Z{t) (8) ^ t )=[l- (1 



q)Pi(t)r 2 



Zi(t), 



(14) 



if \-{l-q)l3{t)r x > 0, and p(r, t) = if l-(l-q)/3(i)r A < 
(cutoff condition) . By replacing this generalized gaus- 
sian function in Eq.(^), we verify that A = 2+9, q = 2 — v, 



with (3x(t) - t-2/[2+d(i-g)] and Zl (t) ~ t d/[2+d(i- q) ]_ Qn 
the other hand, for v = 1, the point-source solution for 
the O'Shaugnessy-Procactia equationQ emerges, 



(3{t) = /3 Q [l + At]- x ^ x+d ^\ 



and 



Z(t) = Z [l+M] 



d/[\+d(l-q)] 



(9) 



(10) 



p(r,t) 



Z2(t), 



(15) 



where j3 2 ~ t' 1 and Z 2 (t) - t d / ( - 2+e \ Moreover, Eq. (§) 
recovers the time-translated Gaussian solution for the Eq. 
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(0), 

P(r,t) 
In fact, Eqs. (0 



(16) 



[4TTD(t Q + t)] d /2 
@ and ([h]) can be respectively written as 

Pit) = i(t +tr A/[A+d(i - 9)1 (17) 



and 



z(t) = B(t + t) d/[A+ < i(1 - <?)1 , 



(18) 



with 1 = /3 A- A /[*+<i(i-<?)] i s = ZoA^/^C 1 -*)], and 
t = I /A. Eq. (|l|) is obtained from Eq. (|) with /3{t) 
and Z(i) given above in the limit v = 1 and A = 2, and by 
using the normalization f2d J p(r,t)r d ~ 1 dr = pa, where 
fid = 27r d / 2 /r(<i/2) is the solid angle in a d-dimcnsional 
space. 

To conclude this review, we notice that our point-like 



solution can be connected with the Tsallis entropy) 20 



III. NONLINEAR FOKKER-PLANCK 
EQUATION WITH TIME-DEPENDENT SOURCE 
TERM 

Nonlinear diffusion with absorption arises in many ar- 
eas of science such as engineering, biophysics, solid state 
physics and reactor, among others Ell E2, 23, EJ. This 



section is dedicated to the analysis of the generalized 
Fokker-Planck equation presented in the introduction 
with a source term. This kind of equation is expected to 
provide a basis for phenomenological applications, where 
a large class of nonlinear anomalous diffusion with source 
plays an important role. 

Let us consider the nonlinear Fokker-Planck equation 
(|l|) . The source term in this equation can be removed by 
an appropriate change in the solution: 



a{t')dt' I /3(r,t). 



p(r, t) = exp 

In this way, /5(r, t) obeys the equation 
dp{r,t) 



(19) 



01 



V ■ {KV[p(r,tW} - V ■ \pF(r,t)p(r,t)], 



(20) 



where K(t) = if exp (1 - v) f* a(t')dt' 



Thus, Eq. 



( |20| ) has the same form of the corresponding one without 
the source term, but with a time-dependent diffusion co- 
efficient K. Note that this time-dependence is induced 
by the nonlinear term p v , disappearing when v = 1. 

As emphasized in the introduction, below Eq. (0), the 
replacement of the operator V • (KV) by A unable us 
also to analyze situations with fractal dimension in the 
sense of O'Shaughnessy-Procaccia. Thus, to investigate 
a radial nonlinear anomalous diffusion equation with a 
source term and non-integer dimension d, we consider 



dp(r,t) 

at 



= DA[p(r,t)Y -a(t)p(r,t) 



(21) 



instead of Eq. ([!]). In the interests of brevity we omit- 
ted the drift term. The case with an external force is 
discussed in Sec. V. To eliminate the source term, we 
proceed analogously as in the previous paragraph. There- 
fore, we obtain 



dpjr, t) 
dt 



DA[p(r,t)Y 



with D(t) = Dcxp (1 - v) f*a(t')dt' 
fine the time variable introducing an effective time as 



(22) 



Now, we rede- 



r(t) = / D(t')dt', 
Jo 



(23) 



so Eq. (t22j) can be written as 



dpjryr) 
8t 



(24) 



The solution of this equation was presented and discussed 
in Sec. II. Thus, it is given by Eqs. (|), (|) and @ with 
D = 1 and t replaced by r. 



Another possible generalized Fokker-Planck is 



d[p(r,t)]i 
dt 



= V • {KV[p(r, t)]"} - V • { M P(r, t)[p(r, *)]"} - a(t)[p(r, t)]\ 

I 



(25) 



This equation, with K and a being constants, was dis- cussed in Ref. |25[. Note also that Eq. (25) reduces 
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to Eq. (P when we replace \p(r, t)] 7 * by p(r,t). In the 
same direction, this kind of generalization can be incor- 
porated in Eq. (|2l|). In this paper, we will not present 
details about these possibilities because their solutions 
can be directly obtained from Eqs. ([!]) or (^l|) by using 
the identification p — p ri . 



IV. SOME KINDS OF SOURCES 

In order to have an enhanced understanding of the rel- 
evance of the source terms in the anomalous diffusion 
equation, we apply the previous general procedure to 
study some specific cases. In this context, the simplest 
situation to be investigated is when a(t) = ao, where «o 
is a constant. We also investigate a more generic source 
term, a(t) = a n t n . 



A. Constant source 

Now, taking into account Eq. (pl|) , we analyze the case 
a(t) = a . Thus, the effective time, Eq. becomes 



r(t) = D 



exp[(l — v)aot] — 1 
(1 - v)ao 



(26) 



This effective time exhibits three different behaviors. For 
(1 — v)aa < 0, the effective time goes exponentially to 
Too = D/[(v— l)ao], i.e., a stabilization occurs in the dif- 
fusive process. For (1 — v)a.Q > 0, there is an exponential 
increase leading to r(t) w Dexp[(l — v)ocot]/[{l — u)ao], 
when t is large. It is interesting to remind that this kind 
of behavior arises due to the nonlinearity of the diffusion 
equation. In this direction, we have r = Dt when v = 1, 
which corresponds to q — 1 as long as q = 2 — v. 

By using Eqs. (0) and (|l(]) with t replaced by r, we 
obtain 



0{r(t)) =0a[l + A 



and 



exp[(l — u)aot] — 1 
(1 - v)a 



Z(r(t)) = Z [1 + A 



exp[(l — u)aot] — 1 
(1 - v)a 



-\/[X+d(l-q)] 



(27) 



d/[\+d(l-q)} 



(28) 



Thus, from Eqs.(||) and (|9|) we verify that 

p(r, t) = [I - (1 - q)/3(r(t))r x ] 1 ^ / Z(r(t)) 

(29) 

if 1 - (1 - q)P{r{t))r x > 0, and p(r,t) = if 1 — (1 — 
q)P(T(t))r x < 0. 



In this solution, when (1 — v)a$ < 0, we have, for large 
time, 



p(r, t) « [1 - (1 - q)^} 1/{1 ~ q) / Z a 

where 

(3 X = /3b {1 + A /[{y l)a ]}- A/[A+d(1 - 9)1 

and 

Z oo = Z {\ + A/[{v-l)a ]} d '^ 1 -^. 



(30) 



(31) 



(32) 



This means that p(r, t) decays (grows) exponentially in 
the time when ao > (ao < 0), and is spatially limited 
for q < 1 as long as p(r,t) = if r A > 1/[(1 — q)(3 00 ]. In 
contrast, the spatial asymptotic power law behavior, 



P(r,t) 



.A/(i-«) 



(33) 



arises for sufficiently large t and (1 — v)o.q > 0. On the 
other hand, when v = 1, we have 



p(r, t) - r d/x cxp[-a t - r x /{DX 2 t)} 



(34) 



for large t. 

Let us focus attention now on the mean square dis- 
placement. By using the effective time in Eq. (|l3|), we 
obtain 



(r 2 ) = C 2 P~ 



2/A 



1 + A 



,(l—u)aat _ i 

(1 - v)ao 



2/[X+d(l-q)] 



(35) 



The long time behavior for (r 2 ) leads to an exponential 
growth when (1 — v)ao > 0. For (1 — v)a.Q < 0, this 
expression converges to an asymptotic value. For v = 1, 
we have (r 2 ) ~ t 2 / x , i.e., when v = 1, the mean square 
displacement is not affected by the constant source term. 
These results show that the presence of a source term 
plays a relevant role on the diffusive process when v ^ 1. 



B. Source with a general power term 

Now we analyze the more general case corresponding to 
the source a n (t) = a n t n , where a n and n are constants. 
The related effective time, T n (t), is given by 



r„(t) 



D exp 



(1 v)ctn /,/n+l __ i n + x \ 

n+1 ^ 1 



dt', (36) 



where to is a cutoff only used to avoid divergence when 
n < —1 and (1 — v)a n /{n + 1) > 0. When n > —1, we 
get 



Here T(c) and T(c, x) are the gamma and the incomplete 
gamma functions, respectively. By using the asymptotic 
behavior for T(c,x), we obtain r n ~ t~ n exp[— a n (y — 
l)i 1+n /(l + n)] for n > — 1 and large t. For n — —1, we 
haver_i(£) ~ fpi-x{i-u)+i ^ and for n < —1, we verify that 
i~ n {t) ~ t- Taking into account the asymptotic behavior 
of the mean square displacement, (r 2 ) ~ r 2/[A+<i(i 
we identify two different behaviors when n > — 1: if 
a n (y — 1) > 0, the mean square displacement reaches a 
limit value, and it grows exponentially if a n (y — 1) < 0. 
For n = -1, its value goes with f2[i-(i-g)a_ 1 ]/[A+d(i-g)] > 
Finally, for n < -1, we obtain (r 2 ) - f 2 /[A+d(i-«)J. 



where we set p, equal to unity. In this section we are 
going to consider a specific case of a particle moving in a 
harmonic spherical potential V(r) = kr 2 /2. Thus, r = 
is a stable equilibrium position for positive k, and r = 
is an unstable equilibrium position for negative k. Eq. 
(|38|), with this potential, corresponds to the Uhlenbeck- 
Ornstein process p6[ in the particular case of 9 = and 
q=l{v = 2-q). 

By using, again, the ansatz (^) into Eq. (pq), we verify 
that Z(£) and /3(i) obey the equations: 



EXTERNAL FORCES 



1 dZ 
~Z~dT 



= DduXpZ 1 *- 1 - kd 



(39) 



Now we draw our attention in order to analyze the 
generalized radial Fokker-Planck equation with a drift 
term. We take Eq. (||) with a radial external force, i.e., 



dp(r,t) 
dt 



= DA[p(r,t)r 



1 d 

r d—l Q r 



„d-l 



F(r)p(r,t)] , 



(38) 



and 



ld0 
]3~dt 



-DvX 2 ^- 1 + kX. 



(40) 



The solution for these two nonlinear coupled equations is 



Z(t) = Z Q 



+ 1 (DuXp^- 1 -k)(l- e -MA+««i- 9 )]t) 



d/[A+d(l-g)] 



(41) 



and 



/3(i) = /?o 



+ i (i^^o^- 1 - k) (l - e -MA+d(i- 9 )]^ 



-\/[\+d(l-q)] 



(42) 



whose initial conditions are Z(0) = Zq and j3(0) = (3q. 

Since A + d(l — g) > 0, for positive fc, we have a sta- 
tionary solution and (r 2 ) goes to a constant value. When 
k < 0, there is no stationary state and the (r 2 ) increases 
exponentially for large t. Our solution with 9 = cor- 
responds to one obtained by Tsallis and Bukman || for 
the one-dimensional case when (x) =0. We also remark 
that O'Shaugnessy-Procaccia equation, Eq. (||), when 
an external force is incorporated, is a special case of our 
results if we take v = 1. 

Before concluding this section, we note that in the pres- 
ence of the source term, — a(t)p, a similar procedure as 
presented in Sec. Ill can be employed, leading to Eq. 
@ with D replaced by D(t). Thus, Eqs. @ and @ 



must be solved with D(t) instead D to obtain the gener- 
alized gaussian solution. 



VI. STATIONARY SOLUTIONS FOR THE 
NONLINEAR FOKKER-PLANCK EQUATION 

We can extend our study about Fokkcr-Planck-like 
equation by investigating stationary solutions. The usual 
Fokker-Planck equation has a stationary solution that 
corresponds to the canonical Boltzmann-Gibbs distribu- 
tion, P oc e~@ v . On the other hand, the corresponding 
stationary solution of the nonlinear equation (0) is [[16] 
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VII. SUMMARY 



p(r)=p [l-{l-q)(3V{v)f^ 



(43) 



if l-(l-g)/3F(r) > 0, and p(r) = if 1 - (1 - g)/?V(r) < 
0. Here = pp 1 ^ 1 /[(2~ q)D], pa is a positive integration 
constant, and F(r) = — W(r). Notice that j3 > when 
V(r) is a confining potential. 

When we consider the unified nonlinear equation 



dp(r,t) 
dt 



V-{i^(r)V[p(r,t)r}-V 



[ /1 (r)F(r)p(r,t)] ) 

(44) 



its stationary solution that generalizes Eq. (|43|) is 

1/(1-9) 



P(r) = Po 



i- 1 ^)aT 1 i ;//(r) 



(45) 



if 1 - (l-g)pg - V e ///(2 - 9) > 0, and p(r) = if 
1 — (1 — q) pl~ 1 V e ff/(2 — q) < 0, where we introduced 
the effective potencial V e //(r) = J ||^VV(r) • dr. Note 
that [p(r) / K (r)]WV (r) must be written as a gradient so 
it is possible to obtain an effective potential in the N- 
dimcnsional case. In particular, for a 3-dimcnsional ex- 
ample, this implies that V[/x(r) / 'K(r)] x VV(r) = 0. This 
condition can be easily accomplished if we take p = p(r), 
K = K(r), and V = V(r). Thus, Eq.@ can be written 



p(r)=p [l-P'(l-q)V eff (r)] 1 /^\ 



(46) 



with 0> = pl 1 /(2 q) and V eff (r) = f $±«ZM dr. 
Following the development presented in the previous sec- 
tions, we stress that the stationary solution (|4^) remains 
true for an arbitrary non-integer dimension. For the one- 
dimensional case, the solution ( |45| ) was obtained in Ref. 

s 



In this paper we have obtained exact spherically sym- 
metric solutions for a nonlinear Fokker-Planck equa- 
tion with radial-dependent diffusion coefficient, time- 
dependent source, and spherical external force. We re- 
mark that the time-dependent solutions are generaliza- 
tions of the gaussian ones. 



Without both external force and source term, there is 
a competition between the fractal aspect and the nonlin- 
ear behavior (dictated by the parameter v), leading to 
superdiffusive, 'normal' and subdiffusive processes. On 
the other hand, in the presence of a source term, the 
mean square displacement is only affected in the nonlin- 
ear case, i.e., the source term modifies the superdiffusive, 
subdiffusive or 'normal' characters of the process when 
v ^ 1. In particular, this behavior is notorious when a 
power law time-dependent source is considered. 

The diffusive process with a external force F(r) = —k r 
leads to different behaviors to the time dependence of dis- 
tribution function. For positive k, we have a stationary 
solution and the mean square displacement goes to a con- 
stant value. When k < there is no stationary state and 
the mean square displacement increases asymptotically 
in an exponential way. 



In the case of the stationary solutions for our Fokker- 
Planck-like equation, we obtained a generalization of the 
canonical Boltzmann-Gibbs distribution. This general- 
ized solution only exists if we are able to obtain an effec- 
tive potential that includes the true potential, the drift 
and the diffusion coefficients. 
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